A quantum cellular automaton for universal quantum computation 
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I describe a quantum cellular automaton capable of performing universal quantum computation. 
The automaton has an elementary transition function that acts on Margolus cells of 2 x 2 qubits, 
and both the "quantum input" and the program are encoded in the initial state of the system. 
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I. INTRODUCTION 

A number of physical systems that are considered for 
the realization of a universal quantum computer, such as 
optical lattices or arrays of micro- lenses 0], possess a 
translation symmetry in the arrangement of qubits and 
their mutual interaction. Quantum cellular automata 
(QCA) represent a suitable framework to explore the 
computational power of such physical systems, because 
they respect this symmetry. A priori, translation invari- 
ance may be regarded as a severe limitation since it con- 
strains the degree of control that can be exerted to the 
quantum system. However, it has been demonstrated 
that one-dimensional QCA can efficiently simulate any 
quantum Turing machine 

Further it has been shown that there exists a univer- 
sal QCA which can simulate any other automaton with 
linear slowdown 0|, and that every reversible QCA can 
be represented in a generalized Margolus partitioning 
scheme 0. Proposals with an emphasis on experimen- 
tal viability have outlined how generic physical systems 
can be used as quantum computation devices if equipped 
with a minimal amount of external control. Among the 
described mechanisms arc global control via sequences 
of resonant light pulses Q or modulation of a coupling 
constant 0, ^iid individual control over one of the 
elementary cells 0- 

At this point, one may abandon all algorithm-specific 
control during the process of computation and ask "How 
intricate do quantum cellular automata have to be such 
that they can perform useful tasks in quantum informa- 
tion processing?". A quick answer may be "Simple, by 
definition." However, when QCAs are tuned for algo- 
rithmic application, it may occur that — while the simple 
composition is retained — the elementary cells and neigh- 
borhood schemes become complicated. An interesting 
facet of the answer to the above question has been pro- 
vided in |lOl |. where a very simple QCA for quantum 
data transmission has been devised (also see Q)- Mo- 
tivated by a recent result where universal compu- 
tation via autonomous evolution of a 10-local Hamilto- 
nian is described, here 1 consider quantum computation 
in the cellular automaton scenario. 1 explicitly construct 
a computationally universal two-dimensional QCA whose 
transition rule is based on a four-qubit unitary. 



II. CONSTRUCTION OF THE AUTOMATON 

Consider a two dimensional lattice of size 2s x 2r with 
periodic boundary conditions, i.e. a torus. Each lattice 
site carries a qubit. The transition rule for the QCA is 
described in terms of a Margolus partitioning [ij. The 
lattice is partitioned into cells of size 2x2, and there is 
a separate partitioning for the time t being even or odd, 
respectively. One may choose a coordinate system on 
the torus with axes parallel to what were the boundaries 
before identification. For t even, the qubits in the upper 
left corner of each cell have both coordinates even, and 
for t odd they have both coordinates odd. Thus, a cell 
in step t overlaps with four cells of step t — 1 . 

The transition of the QCA from time tiot+ \ proceeds 
by simultaneously applying a unitary transition function 
T to each cell. For a suitable choice of the 4-qubit unitary 
T one can perform universal quantum computation with 
the described QCA. Specifically, a quantum logic network 
of local and next-neighbor unitary gates with width 2s 
and depth r can be simulated. 

With a labeling of particles as illustrated in Fig. 
the following elementary transition function is chosen: 
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Therein, S(a, b) denotes a SWAP-gate between qubit a 
and b, Hi is a Hadamard transformation on qubit 1, and 
Zc denotes a Pauli phase flip operator applied to qubit c. 

If 1^)34 is a state in the computational basis, it effec- 
tively stores two classical bits, p*-^-* and p^^-* . Then, the 
transition function r amounts to a classically controlled 
unitary operation U{p) applied to |£')i2- The bit p^^-* 
triggers a A{Z) gate applied to 1-0)12, and p^^^ a 7r/4- 
phase gate exp(— Itt/S Zi) applied to qubit 1 of |-D)i2- In 
this way, \p) encodes an elementary step of a program, 
carried out on the two-qubit "data" |£')i2- The SWAP- 
gates allow the quantum data and the program to pass 
by another, such that |-D)i2 can interact with subsequent 
program registers. 

The sets of qubits on the torus with the same first 
coordinate x, < a: < 2r — 1, are called columns. At 
time t = all even columns contain data registers \Di), 
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FIG. 1: Transition of the QCA from t to a: Margolus cell 
of 2 X 2 qubits, with the part |-D)i2 of a data- and |p)34 of the 
program column, b: Before and after the first transition. The 
program columns move left and the data columns move right. 
The dashed lines indicate the partitioning into Margolus cells. 



and all odd columns contain program registers \pj) 
initial state of the automaton is 
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FIG. 2: The quantum logic networlt simulated by the QCA. 
The gates in each shaded box result from one application of 
the elementary transition function. The arrows denote clas- 
sically controlled gates which are triggered by the program 
registers. 



Let Ti,j denote the elementary transition function 
applied to the cell Therein, i is the column co- 

ordinate of the upper left qubit in the cell, and j the 
respective coordinate within the column. Then, the uni- 
tary transition function Tj, acting upon two consecutive 
columns i, i + 1, is 



1^(0)) =(g)|A(0))2.b.+ l)2.+l• 



(2) 



Here and in the following the labels inside the kets specify 
the state and the ones outside specify the location of the 
support within the lattice i.e., the column. For example, 
|^i(0))2 is data register No. 1 at time t = 0, located 
on column 2. Of all data registers only |Do) is used, the 
others arc auxiliary. When the QCA starts to run, the 
data registers move right (counter-clockwise, as seen from 
top) and the program registers move left (clockwise), by 
one column in each time step. When passing the data 
registers, the program registers \p) control unitary trans- 
formations U{p) acting upon the data registers \D). In 
this way, a program specified by the data pi, ■■■,Pr-, with 
Pi encoding the first and pr the last step, is carried out on 
the quantum data register \Do). The program that is car- 
ried out corresponds to a quantum logic network of local 
and next-neighbor gates in a particular arrangement; see 
Fig. 121 Such networks are sufficient for universal quan- 
tum computation, as is discussed in detail further below. 
The same program steps that are applied to the register 
I Do) E^rc also carried out on the auxiliary data registers 
\Di), 1 < i < r ~ I, but in scrambled order. Therefore, 
these registers are not used. 

In course of computation, both data and program 
travel across half the torus. When the automaton has 
run for r time steps, the computation is finished and the 
register \Do{r))r can be read out from column r, via local 
measurements. 



T, 



(3) 



Therein, [2j + i]2s is a shorthand for 2j + i mod 2s that 
will be used throughout the remainder of the paper. 
If \p) is a state in the computational basis, then 



(|i?).®|p).+i) = |p).®(f/(p)|i?)),+i 



(4) 



Therein, U (p) is a unitary transformation chosen by p 
containing a Hadamard- and possibly a A(Z)- and a 7r/4- 
phase gate, in accordance with . 

The global transition function T : \tp{t)) — > \ip{t+l)) 
is, for even t given by Tc = <S>lZo ^2»' ^^'^ 0*^^ ^ ^y 
To ~ ^1=0 -^2i+i- In both cases it can be written in the 
form 



T 



[2i + t]2r- 



(5) 



i=0 



Now, the state |V'(^)) of^ the QCA at time t is 



i=0 



with the data register i at time t given by 



|A(t)) 



|A(0)). 



(7) 
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The unitarics C/(p[i+fe]^) arc ordered in ascending order 
with k, i.e. U(jj[i^i^^) acts first. For i ~ and t = r, in 
particular, one finds that 



1^0 (r)) 



n u^p^) 



.fe=i 



l^o(O)) 



(8) 



is the output quantum register, with the unitaries 
U{pi)...U{pr) applied in the correct order to the quantum 
register in its input state, |Do(0))- 

Before proving ©, let us recover therein some fea- 
tures of the QCA that were stated before. It is easy 
to see that — apart from being moved — the program reg- 
isters remain unchanged throughout the evolution, and 
that there is no entanglement across columns. The data 
registers indeed move right, and the program registers left 
{withi+t+l = i', |p[j+t+i]J[2i+t+i]2, = b[j']J[2,'-t-i]J. 
After r time steps, the output quantum register \Do{r)) 
can be read out from column r. 

Eq. lO is proved by induction. First note that for 
t = Q it reduces to Further, 

r-1 

T\lp{t)) = (^T^2^+t]2r\D^{t))l2^+t]2,^Pl^+t+l]^)l2^+t+l]2r■ 
i=0 
r-1 

\P[^+t+l],^[2^+t]2r■ 

r-1 

= (8)|A(<+l))[2,+ («+i)],,® 



b[i+(t+i)+i].)[2i+(t+i)+i]2 
H'{t + i)). 



(9) 



Here, the first line follows by the definitions of T and 
© and lO, and the second follows by ©. The 
third line follows by {Tj) and, for the |p)-part, the substi- 
tution i ^ i + I under which the product is invariant. □ 

Finally, it is shown that the quantum logic network 
simulated by the described QCA is indeed universal, as 
stated. The CNOT-, the Hadamard- and the 7r/4-phasc 
gate exp(— i7r/8Z) form a universal gate set For 
the described QCA, one can independently apply the 
Hadamard-, the 7r/4-phase gate and the identity to the 
simulated logical qubits, by choosing the following se- 
quences of p*^'^^ -bits 



000 000 000, for U 
^(3) - J 0101101101, for U 
1000 000 000, for U 



I, 
H, 

C/.[7r/4]. 



(10) 



with all p^^-'-bits zero. In the way they are constructed 
here (no claim that this is any close to optimal) , the one- 
qubit operations from the universal set require 10 succes- 
sive applications of p''^) -controlled gates and thus 20 time 
steps of the QCA. A next-neighbor CNOT-gate that acts 
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FIG. 3: Boundary specifications for the described QCA on 
a planar sheet. The 2 x 2-cells displayed in white are acted 
upon by r, as usual. To the hatched 1 x 2-cells a SWAP- 
gate is applied, and to the cross-hatched cells the identity 
operation. The column in gray underlay represents a data 
register moving right. 



within this cycle can be constructed from a A{Z)- and 
two Hadamard-gates. The long-distance CNOT-gates 
may then be constructed with the help of next-neighbor 
SWAP-gates, which themselves consist of three next- 
neighbor CNOT-gates. This completes our construction 
of a QCA capable of performing universal quantum com- 
putation. 

Two remarks: 1) The described QCA may, with some 
right, be called a deterministic programmable quantum 
gate array, but this notion is already in use for a construc- 
tion that has been proven not to exist 0|. Our QCA is 
consistent with this result. The program information is 
classical and all program states {S)i=i \Pi)2i-i are orthog- 
onal, as required in [l3 |. Further, from the viewpoint 
of temporal complexity, the described QCA is — within a 
constant — as efficient as a quantum logic network with 
local and next-neighbor gates. 

2) That the described QCA lives on a torus simplifies 
the discussion, but is not essential. A planar sheet of 
size n X 2r is sufficient for simulation of the discussed 
networks of r time steps and n qubits. Fig. |2| specifies 
the operations on the boundary that differ from r. When 
this modified QCA is run, in the bulk the data registers 
still move right and the program registers left. On the left 
and right boundary, however, the registers are reflected. 
As a consequence, on the left side of the lattice, reflected 
program registers are acted upon by left-moving program 
registers as if they were data. More severely, on the right 
side of the cluster, reflected data registers act upon right- 
moving data registers as program. Therefore, the state 
of the QCA is no longer a tensor product of the column 
states, but instead an entangled state supported by many 
columns grows from the right, by one column in each time 
step. However, none of that has an impact on the data 
register \Do), which — as before — can be read out from 
column r after r transitions of the automaton. 
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III. CONCLUSION 

I have described a quantum cellular automaton for uni- 
versal quantum computation. The transition function 
from one time-step to the next is generated by a four- 
qubit unitary transformation acting on Margolus cells of 
size 2x2. The program is encoded in the initial state of 
the system, and the automaton is left to its autonomous 
evolution from initialization to readout. 

Coming back to our initially posed question, it is found 



that QCAs performing complex tasks in quantum infor- 
mation processing can indeed be constructed for compact 
cells and neighborhood schemes. 
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